Quantum correction to the Kubo formula in closed mesoscopic systems 
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We study the energy dissipation rate in a mesoscopic system described by the parametrically- 
driven random-matrix Hamiltonian H[ip(t)] for the case of linear bias tp = vt. Evolution of the field 
ip(t) causes interlevel transitions leading to energy pumping, and also smears the discrete spectrum 
of the Hamiltonian. For sufficiently fast perturbation this smearing exceeds the mean level spacing 
and the dissipation rate is given by the Kubo formula. We calculate the quantum correction to the 
Kubo result that reveals the original discreteness of the energy spectrum. The first correction to 
the system viscosity scales tx v~ 2 ^ 3 in the orthogonal case and vanishes in the unitary case. 
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The Kubo formula § is one of the cornerstones of 
modern condensed matter physics. It is a standard tool 
for calculating various linear response functions, with 
conductivity as a prototypical example. Based on the 
lowest-order perturbation theory for a continuous spec- 
trum, Kubo approach connects the kinetic response of a 
system with respect to some external field <p(t) to the 
equilibrium correlator of generalized currents. 

Application of the Kubo formula essentially relies on 
the assumption of a continuous spectrum S. (For the 
case of a discrete spectrum it gives either zero, when the 
frequency of the external field ip(t) is out of resonance 
with any pair of energy levels, or infinity otherwise, in- 
dicating breakdown of the linear response theory.) The 
widespread usage of the Kubo formula for macroscopic 
objects is justified by the smallness of the mean level 
spacing A oc H d / L d (d is the dimensionality of space and 
L is the system size) which is usually smaller than the 
inelastic width T[ n of energy levels. In an electron sys- 
tem, the inelastic smearing is due to interaction as well 
as escape to reservoirs operative for open systems. 

For closed microscopic systems (e.g., quantum dots), 
the condition of continuous spectrum is violated at suf- 
ficiently low temperatures when the interaction-induced 
smearing Ti n becomes smaller than A. In this case in- 
dividual discrete levels are well resolved and the Kubo 
formula may fail to describe the dissipation in the sys- 
tem. Such a situation had recently been discussed in the 
context of vortex dynamics in layered superconductors at 
low temperatures ||[ [| |[ ^) . 

In a closed system with Fi n <C A, the mechanism of 
dissipation depends on the rate v = d<p/dt of variation 
of the external field <p(t). For sufficiently slow pertur- 
bations with v <C vk (the velocity vk depends on the 
sensitivity of the spectrum to the change of ip and will 
be defined below), the system adiabatically follows the 
spectrum of its instantaneous Hamiltonian H[ip(t)], and 
the dissipation is due to rare Landau-Zener transitions 
between individual levels H. Fast perturbations with 
v ^> vk cause multiple transitions between energy levels 
thereby transforming the discrete spectrum of the sta- 
tionary Hamiltonian into a featureless quasi-continuous 



spectrum, where the dissipation can be obtained with the 
help of the Kubo formula. 

These two regimes had been studied by Wilkinson 
who calculated the dissipation in the adiabatic (v <C vk) 
and Kubo (v 3> vk) regimes assuming the Hamiltonian 
from the Wigner-Dyson random- matrix ensembles || . In 
the Kubo regime, the energy dissipation rate is given by 
the linear-response formula: 



W Kuho = ^C(0)V 2 , 



(1) 



where (3 — 1 for the orthogonal (GOE) and 2 for the uni- 
tary (GUE) ensembles, and C(0) = ((dEi/dcp) 2 ) /A 2 is 
the variance of the level velocity normalized by the mean 
level spacing A. (The system of units with h = 1 is used 
throughout the Letter.) Equation ([!]) suggests that C(0) 
has the meaning of the generalized conductance [|l0). It 
determines the critical velocity vk ~ A/ y/ C(0) separat- 
ing the adiabatic and Kubo regimes of dissipation. In the 
adiabatic limit, the dissipation rate nontrivially depends 
on the symmetry of the Hamiltonian B: 



(2) 



where c x = (7r/4)r(|)[C(0)] 3 / 4 A 1 /2 and c 2 = ttC(O). 
Hence the dissipation is superohmic for the GOE, while 
for the GUE it remains ohmic, exactly coinciding with 
WKubo despite a very different mechanism of dissipation. 
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FIG. 1: Position of the adiabatic and Kubo regimes as a func- 
tion of the inelastic width Ti n and velocity v. 
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A schematic diagram indicating the regions of the adi- 
abatic and Kubo regimes as a function of the inelastic 
width Tin and velocity v of external perturbation is shown 
in Fig. [l| Counterintuitively, the linear-response Kubo 
formula does not describe the low-velocity behavior of 
closed systems at r in <C A! 

In this Letter we address the question of how does 
the discreteness of energy levels of a stationary (v = 0) 
Hamiltonian manifest itself in the Kubo regime (v 3> vk) 
when the levels are smeared into a quasi-continuous spec- 
trum. We show that the relative correction to the high- 
velocity asymptotics (Q) can be regularly expanded in 
integer powers of (vk/v) 2 ^ 3 . For the orthogonal symme- 
try, the first (one-loop) quantum correction to Wxubo is 
given by 
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(3) 



where the omitted terms have the order of 0[(vk /v) 4 / 3 ], 
and the crossover velocity is defined as % = 
(2 1 / 2 /7r 2 )A/ v / C(0). Thus, the remaining correlations 
in the quasi-continuous spectrum enhance dissipation at 
t; > njf, with a gradual crossover to the superohmic 
regime (^) at u C ji^. In the unitary case, the first and 
the second (two-loop) quantum correction to the Kubo 
result ([!]) vanish, making it tempting to conjecture that 
the identity Wi = Hi<ubo holds for all velocities. 

Quantum corrections to the quasiclassical properties 
of disordered systems had been the subject of intense 
studies in the last decades [jll], |l2). In treating these 
phenomena, the nonlinear er-model was proven to be the 
most powerful tool both in the perturbative and nonper- 
turbative regimes (where it often is the only possible ap- 
proach) . Three versions of the a- model based on the su- 
persymmetry [|12| , replica , and Keldysh |l4], [l5) tech- 
niques had been proposed for noninteracting systems. 

The problem of energy pumping by the parametrically- 
driven Hamiltonian H[<p(t)] belongs to the class of non- 
equilibrium problems, that dictates the choice of the 
Keldysh formalism as a solution tool. We will derive the 
Keldysh cr-model for the parametrically-driven random- 
matrix Hamiltonian and show that its saddle-point solu- 
tion yields the kinetic equation for the distribution func- 
tion, reproducing the Kubo result ([!]) . Fluctuations near 
this saddle point are responsible for the quantum correc- 
tion to the Kubo formula leading to Eq. (||). 

We consider a time-dependent matrix Hamiltonian 



H[ip(t)] = Hi cos kip(t) + i?2 sin k(p(t), 



(4) 



where Hi and Hi are the N x N matrices from the 
same GOE (H T = H) distributed with the probabil- 
ity density P(V) cx exp[-(7r 2 /2JVA 2 ) trV 2 ], and k = 
(7ro-/A)(2/A r ) 1 / 2 . The density of states for an instant 
Hamiltonian is given by the Wigner semicircle: p{E) = 
A _1 [l - tt 2 E 2 /4N 2 A 2 ] 1 / 2 , with A being the mean level 



spacing at the center of the band. The dispersion of the 
matrix elements 



dm 



dtp 



= a 2 (l + < %) 



(5) 



is independent of the matrix size N. The generalized 
conductance is then C(0) = 2a 2 / A 2 . The states of the 
Hamiltonian are filled by N/2 noninteracting fermions. 

The time evolution governed by ip(t) will change the 
state of the system and eventually pump some energy 
into it. The energy of the system will grow unless the 
inelastic interaction with the thermal-bath is taken into 
account. This interaction will establish a nonequilibrium 
steady state. Remarkably, however, that the energy dis- 
sipation rate is independent of the resulting distribution 
and hence can be calculated as a mean growth rate of the 
total system energy of noninteracting fermions j^] . 

Within the Keldysh formalism, the system is described 
by the action 



sm = -i f 

J— t 



dt W{t) 



(6) 



where ^(t) is a Grassmanian 47V- vector field acting in 
the direct product of the index space of the Hamiltonian, 
Keldysh (K) space and particle-hole (PH) space. The 
Pauli matrices in the K and PH spaces are denoted by 
<7i and Ti, respectively. The Keldysh forward-backward 
structure is a necessary ingredient of the dynamic formal- 
ism ensuring the correct normalization Z = j D^e^ s — 
1 in the absence of quantum sources, whereas the PH 
grading is introduced in order to handle the orthogonal 
symmetry of the Hamiltonian Jl^| . 

Derivation of the a-model, which is a low-energy effec- 
tive theory for the action (^), is a standard procedure. 
One has to average Z over the Hamiltonian (^J), decouple 
it by the 4x4 matrix Q u > and integrate over fermions \P. 
Keeping the terms which are finite in the limit N — » oo 
and assuming linear bias <p = vt we arrive at the following 
action (the weight is e~ s ) for the cr-model 16 : 



7" m A ^ 7T^ 3 



dtdt' (t-t' ftri 



t, (7) 



where il 3 = ircr 2 v 2 /A = (tt/2)C(0)w 2 A. The first term in 
Eq. (^) is the standard random-matrix action |jl7| which 
is responsible for the whole spectral statistics. The sec- 
ond term is of kinetic origin; it accounts for interlevel 
transitions of the time-dependent Hamiltonian H[ip(t)]. 
The field theory with the action ([?]) describes both the 
adibatic and Kubo regimes of dissipation on an equal 
footing. It is controlled by the single parameter 
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A 



v 

VK 



2/3 



(8) 



which will be used hereafter as a measure of velocity v. 



In the stationary case [CI = 0), the Keldysh Green 
function Q is diagonal in the energy representation p5[: 



(a) 




(b) 
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(9) 



where F(E) = 1 - 2f(E), and f(E) is the electron distri- 
bution function. The evolution of the distribution func- 
tion at SI ^ is described by the saddle point of the 
action (pj). Varying with respect the constraint Q 2 = 1 
one obtains the saddle-point equation [Q, SS/SQ] = 0. 
Seeking the solution in the form (^) , we obtain the equa- 
tion for the distribution function f t t'- 



d_ 

dt 



d_ 

or 



f tt , = -n 3 ( t 



t'ffw. 



(10) 



The same equation had been obtained in Ref. [T^ start- 
ing with the diagrammatic technique. Performing the 
Wigner transformation f(E 7 t) — J dre lEr ft+ T /2,t~r/2 
we arrive at the kinetic equation 



df(E,t) =n3 a 2 f(E,t) 



Of 



dE 2 



(11) 



This is a diffusion equation in the energy space, with 
the bare "diffusion coefficient" T)q = SI 3 . The rate of 
energy pumping for the system described by the kinetic 
equation (0) is given by 



W = E 



,d f(E,t) dE _ SI 3 

at "a" ~ a 



(12) 



that coincides with the result of the Kubo formula (pj). 

The kinetic equation ( jll| ) is a true saddle-point of the 
action for all velocities v. The answer (|l2|) for the dissi- 
pation rate predicted is valid, however, only in the Kubo 
regime, and is completely wrong in the adiabatic regime. 
The reason is that the saddle-point approximation is jus- 
tified by the large value of the parameter fl/A. In the 
adiabatic regime, the saddle-point approximation fails 
and one has to take all the Q- manifold into account. As a 
result of this procedure, the answer (||) should be repro- 
duced. Note the interchange of steps with respect to the 
Wilkinson's derivation ||]: He first calculates the proba- 
bility of the Landau-Zener transition and then averages 
it over the distribution of avoided crossings. Here we, 
instead, first average over randomness and then extract 
the dissipation rate from the field theory (0). Thus it is 
a challenging problem to demonstrate how the adiabatic 
result (||) should be obtained from the field theory (]?]). 

Quantum correction to the Kubo result ( |l2] ) in the 
limit SI/ A 3> 1 (A/SI is the loop expansion parameter) 
can be obtained in the regular way by expanding over 
Gaussian fluctuations near the saddle point. The matrix 
Q can be written as |jla 



= u- x pu f , 



u F = 



1 F 
-1 



(13) 



FIG. 2: 1-loop corrections to the system energy. The solid 
fines stand for the propagators ( |l5| ) and the dashed line 
denotes for the source field n(t), and the open and solid ver- 
texes come from the terms S K and S, respectively. 
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FIG. 3: 2-loop corrections to the system energy in the unitary 
case. 



with the Hermitian matrix P having an additional sym- 
metry P T = <j\TiPtiG\ imposed by the orthogonal sym- 
metry of the Hamiltonian. We choose the standard ra- 
tional parameterization of the P-matrix, P = <jj,t^{1 + 
V/2)(l - V/2)- 1 , with the unit Jacobian. The matrix V 
is explicitly given by 



V = 





-at 
-&t 




a o 





b* -c 



\ 

-b T 

J 



(14) 



where the inner (outer) grading corresponds to the PH 
(K) space. The elements au> and btf describe cooperon 
and diffuson modes, respectively. Their bare propagators 
read [fill fill M M 



{ a tit 2 a t3 tJ 



2A r/ 

7T 



2A r/ 

<H*13 - ^24 

7T 



)0(ti 3 )e-° (t " /3+ ' 13t "\ 
(15) 

)9(t 13 )e- Q3t ^, (16) 



ti — tj. In the stationary case (SI = 0), 



(Kt 2 bt 3 t4) 
where tjj 

Eqs. ( ^5[ ) and (16) describe the standard cooperon and 
diffuson in the time domain. The exponential decay of 
the correlators at the time scale SI -1 is a manifestation of 
dephasing by the time-dependent perturbation [|l8| [l9| . 

The system energy (E) (apart from an additive con- 
stant) can then be obtained as a functional derivative 
with respect to the quantum source field n{t): 



(E(t)) = \ §Z[K] 



2 6n{t) 



K=0 



Z[k] = J , 



-S-S K 



DQ, (17) 



where the source action S K — (n/2A) Tr kEo\t^Q. In- 
deed, the presence of the Pauli matrix a\ cuts the off- 
diagonal in the Keldysh space elements of the matrix Q, 
thus projecting onto the upper right (Keldysh) block con- 
taining the distribution function renormalizcd by quan- 
tum fluctuations (the lower left block vanishes due to 
causality) . 
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At the saddle point, Eq. fllTp reproduces the result 
(|l2|). The one-loop diagrams are shown in Fig. ||. 
The diagram (a) coming from the expansion of the 
sourse term S K contains either {a(ti,t 2 )a* (ti,t 3 )) or 
{b(ti 7 t 2 )b*(ti,t3)) which are proportional to 0(0). In the 
Keldysh formalism, the Heaviside ^-function of zero argu- 
ment evaluates to zero, that is related with the causality 
of the theory |L7| . Calculating the diagram (b) we obtain 
for the quantum correction to the dissipation rate: 



5Wi T(i) A 



(18) 



which after rewriting in terms of velocity leads to Eq. (^). 
This result may be interpreted in terms of renormal- 
ization of the bare "diffusion coefficient" in Eq. (JTT]) : 
V — > T>o(l + SWi/WKubo), bearing a natural analogy 
with the weak localization phenomena. In our case, the 
ratio fi/A plays the role of the dimensionless conduc- 
tance as it controls the expansion in terms of the diffusive 
modes. 

In the unitary case, the diagram of Fig. ||(b) vanishes 
indicating the absence of the one-loop quantum correc- 
tion to the Kubo result. In the two-loop approximation, 
only the diagrams shown in Fig. |3| contribute to the dis- 
sipation rate for the GUE. Each of them has the order of 
(A/SI) 2 , but their sum is zero. Thus, for the case of the 
unitary symmetry, the two-loop correction also vanishes. 
Taking into account the coincidence between the low- and 
high- velocity asymptotics (|]) and (||) for (3 — 2, such a 
cancelation is a strong argument in favor of the exact 
identity W% = Wxubo valid for all velocities. At present 
we cannot prove this conjecture but we hope that this 
can be done by an accurate analysis of the cr-model (0). 
We conjecture that the independence of the dissipation 
on the velocity v might be another manifestation of the 
peculiar properties of the unitary group |2^] . It is worth 
mentioning that the direct quantum-mechanical calcula- 
tion of the transition rates at v ~ Vk prior to disorder 
averaging seems completely impossible. Therefore, the 
identity Wi = Wjcubo for the averaged dissipation rate 
would be a highly nontrivial fact. 

The results obtained may be relevant for the descrip- 
tion of heating effect in closed quantum dots whose 
shape is being changed by the low-frequency external 
perturbation. Vortex motion in impure superconduc- 
tors is another field of application, where the conjecture 
W<i = WKubo would indicate the independence of the dis- 
sipative flux-flow conductivity on the velocity of vortex 
motion (|[ g. 

Summarizing, we have developed the Keldysh cr-modcl 
approach to study energy pumping in the parametric- 
ally-driven random-matrix ensembles. We calculated the 
leading quantum correction to the high-velocity dissipa- 
tion, which reveals the original discreteness of the spec- 
trum of the stationary Hamiltonian. This correction 



emerges in the one-loop approximation for the GOE, and 
is absent within the two-loop accuracy for the GUE. 

I am grateful to D. M. Basko, Ya. M. Blanter, 
D. A. Ivanov, M. V. Feigel'man, V. E. Kravtsov, 
A. I. Larkin, and Yu. V. Nazarov for illuminating dis- 
cussions. Financial support from the SCOPES program 
of Switzerland, the Dutch Organization for Fundamen- 
tal Research (NWO), the Russian Foundation for Ba- 
sic Research under grants 01-02-17759 and 02-02-06238, 
the program "Quantum Macrophysics" of the Russian 
Academy of Sciences, the Russian Ministry of Science, 
and the Russian Science Support Foundation is acknowl- 
edged. 



[1 
[2 

[3; 

[4 

[5 

[7: 

[8 
[9 

[10 

[11 

[12. 

[13 



[14 
[15 

[16 



[17 

[is: 

[19 

[20 
[21 
[22 



R. Kubo, Can. J. Phys. 34, 1274 (1956). 

G. D. Mahan, Many-particle physics, 3rd ed. (Kluwer, 

Boston, 2000). 

M. V. Feigel'man and M. A. Skvortsov, Phys. Rev. Lett. 
78, 2640 (1997). 

A. I. Larkin and Yu. N. Ovchinnikov, Phys. Rev. B 57, 
5457 (1998). 

A. A. Koulakov and A. I. Larkin, Phys. Rev. B 60, 14597 
(1999). 

M. A. Skvortsov. D . A. Ivanov, and J. Blatter, preprint 
cond-mat/0207500[ 

L. D. Landau, Phys. Z. Sowjetunion 2, 46 (1932); 
C. Zener, Proc. R. Soc. London, Ser. A 137, 696 (1932). 
M. Wilkinson, J. Phys. A 21, 4021 (1988). 
M. L. Mehta, Random Matrices (Academic Press, Bos- 
ton, 1991). 

B. D. Simons and B. L. Altshuler, Phys. Rev. B 48, 5422 
(1993). 

L. P. Gor'kov, A. I. Larkin, D. E. Khmelnitskii, Pis'ma 
v ZhETF 30, 248 (1979) [JETP Lett. 30, 228 (1979)]. 
K. B. Efetov, Super symmetry in Disorder and Chaos 
(Cambridge University Press, New York, 1997). 
F. J. Wegner, Z. Phys. B 35, 207 (1979); K. B. Efetov, 
A. I. Larkin, and D. E. Khmelnitskii, Sov. Phys. JETP 
52, 568 (1980). 

M. L. Horbach and G. Schon, Ann. Phys. 2, 51 (1993). 
A. Kamenev and A. Andreev, Phys. Rev. B 60, 2218 

(1999) . 

The general form of the kinetic term valid for an arbitrary 
dependence of <p(t) is S kin = (/3tt 2 C(0)/16) / dtdt' [<p{t)~ 
tp{t')] 2 tr Qtt'Qt't- In this form it is valid both for the 
GOE and GUE. 

A. Altland and A. Kamenev, Phys. Rev. Lett. 85, 5615 

(2000) . 

V. I. Yudson, E. Kanzieper, and V. E. Kravtsov, Phys. 
Rev. B 64, 045310 (2001). 

M. G. Vavilov and I. L. Aleiner, Phys. Rev. B 60, R16311 
(1999). 

X.-B. Wang and V. E. Kravtsov, Phys. Rev. B 64, 033313 

(2001) . 

There is also a finite average {&t 1 t 2 &t3t 4 ) but it is irrele- 
vant for the calcu lation of the mean dissipation rate. 
M. R. Zirnbauer, |sond-mat/9903338|. 



